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We present results for inclusive 7r° production in proton-proton and in Au-Au at RHIC energy 
y/s — 200 GeV. We use next-to-leading order perturbative QCD calculation and we include nuclear 
effects such as parton energy loss and nuclear shadowing. We consider the ratio of it distribution in 
Au-Au and p-p collisions for pr > 3 GeV for three cases of parton energy loss: 1) constant parton 
energy loss per parton scattering, = const, 2) Landau-Pomeranchuk-Migdal energy-dependent 
energy loss, e„ ~ V~E n and 3) Bethe-Heitler energy-dependent energy loss, ~ E%. We show 
that recently observed suppression of n° production in Au-Au collisions at RHIC, which is found to 
increase with pt increasing from 3GeV to 8GeV, can be reproduced by e£ = 0.06-E^. We show that 
the ratio of prompt photons to neutral pions produced in Au-Au collisions at RHIC has a strong pr 
dependence approaching one at pr ~ lOGeV. 



Recent results from the Relativistic Heavy-Ion Collider 
(RHIC) showing large suppression of ir° production in 
Au-Au collisions relative to proton-proton collision, in 
the large pr region (3GeV < pr < 8GeV), have created 
new excitement in the field [1]. Since partons produced 
in hard parton-parton collisions propagate through the 
hot and dense medium created in the heavy-ion collision 
and lose their energy, production of large pr pions could 
potentially provide a unique opportunity for studying the 
properties of the hot and dense matter and the possible 
formation of a new phase, the quark-gluon plasma [2-4] . 
Furthermore, parton energy loss is expected to suppress 
the spectrum of the final state pions produced in heavy- 
ion collision as compared to pions produced in hadronic 
collisions [5]. The observed dependence of the sup- 
pression was quite unexpected [1] and posed a great chal- 
lenge for theoretical models which typically predicted a 
rise of this ratio with p t at large transverse momentum 
[5]. 

The aim of this letter is to show that the observed sup- 
pression can be reproduced in an economical way using 
next-to-leading order pQCD calculation of inclusive 7r° 
production [9] augmented with energy dependent energy 
loss and nuclear shadowing. We emphasize that NLO 
pQCD [9] predictions for 7r° production in proton-proton 
collisions at RHIC energy of y/s = 200 GeV was found to 
be in good agreement with PHENIX data [1] for pt > 3 
GeV without a need for introduction of parton intrinsic 
transverse momentum. We show that assuming that 6% 
of parton energy is lost per each parton scattering, one 
can reproduce the observed pt dependence of the sup- 
pression. Here we also present results for prompt photon 
production in proton-proton and in Au-Au collisions at 
RHIC energy of y/s = 200 GeV and discuss the pr region 
in which prompt photons become dominant. 

In perturbative QCD, the inclusive cross section for 
pion production in a hadronic collision is given by: 



E w -^-(\/s,Pn) = J dx a J dx b J dzy^Fi(x a ,Q 2 



F 3 {x b ,Q 2 )D ch {z,Q})E f & ^ X (1 

where Fi(x,Q 2 ) is the i-th parton distribution in a nu- 
cleon, x a and x b are the fractional momenta of incoming 
partons, D c / 7T (z,Q 2 ) is the pion fragmentation function, 
z is the fraction of parton energy carried by the pion 

and d a ^ aX are parton-parton cross sections which in- 
clude leading-order, 0{a 2 ), subprocesses and the next- 
to-leading order, O(a^), subprocesses. 

The parton distribution functions Fi(x,Q 2 ) are mea- 
sured in Deep Inelastic Scattering experiments at HERA 
[6], while fragmentation functions, D c / 7T (z, Q 2 ), that de- 
scribe the transition of the partons into the final-state 
pions are extracted from e + e _ annihilation data from 
PETRA, PEP and LEP [7]. We use MRS99 parameteri- 
zation of nucleon structure functions [8], which includes 
NLO corrections, BKK pion fragmentation functions [7], 
and we set renormalization, factorization and fragmen- 
tation scale to be equal to pt- Varying scale from pt 
to 2pT does not affect the shape of the distribution, but 
rather the overall normalization [9,10]. Our results for in- 
clusive 7T° production in proton-proton collision at RHIC 
energy is shown in Fig. 1. Recent PHENIX data on 7r° 
production in proton-proton collision at y/s — 200GeV is 
found to be in agreement with NLO pQCD predictions 
[1]. Since we restrict ourselves to pt > 3 GeV region, 
we do not include the so called parton intrinsic momenta 
which are expected to be rather small in the high p t re- 
gion [11]. 

To calculate the inclusive cross section for pion pro- 
duction in heavy ion collisions, we will use (1) with the 
distribution and fragmentation functions appropriately 
modified to include nuclear effects such as shadowing and 
energy loss. The modification of the parton distribution, 
known as nuclear shadowing (for a review see [12]), can 
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be written as 

F a/A (x, Q\b t ) = T A (b t ) S a/A (x, Q 2 ) F a/N {x, Q 2 ) 

where T A (bt) is the nuclear thickness function, 
F a / N (x,Q 2 ) is the parton distribution function in a nu- 
cleon and S a / A (x,Q 2 ) is the parton shadowing function. 
We use recent parameterization of the shadowing func- 
tion of Eskola, Kolhinen and Salgado (EKS98), which is 
Q 2 dependent, distinguishes between quarks and gluons 
[13] and is shown to be in very good agreement with the 
NMC data on Q 2 dependence of F 2 Sn /F 2 c [14]. 

We also include the medium induced parton energy 
loss effect. Fast partons produced in parton-parton col- 
lision propagate through the hot and dense medium and 
through scatterings lose part of their energy [15] and then 
fragment into hadrons with a reduced energy. While a 
dynamical study of the parton propagation in a hot and 
dense medium created in a realistic heavy-ion collision 
and the modification of the hadronization is more desir- 
able, we will use a simple phenomenological model [16] 
here to demonstrate sensitivity of pion production to the 
parton energy loss. Given the inelastic scattering mean- 
free-path, A a , the probability for a parton type a to scat- 
ter n times within a distance AL before it escapes the 
system is assumed to be given by the Poisson distribu- 
tion. The modified fragmentation function is given by 
[16], 



JY 



zD a/7V (z,AL 1 Q 2 ) = J2Pa(n) 



(2) 



where z a n = zE%/E%, E a n+l = E° 



Pr/e-m and 



D i is the hadronic fragmentation function which gives 

a/7T ° ° 

the probability that quark or a gluon would fragment 
into a pion. The average number of scatterings within a 
distance AL is (n a ) — AL/X a . We take X a = lfm and 
AL = R A . 

We calculate the invariant cross section for ir° pro- 
duction in heavy-ion collision normalized to the number 
of binary nucleon-nucleon collisions, N co u, where N co a 
can theoretically be determined from nuclear overlapping 
function, i.e. N coll = T AA (b)ag%. The number of N-N 
collisions depends on the centrality that experiment trig- 
gers on. Here we take N co u — 975, which is obtained 
by PHENIX for their central collisions [1]. In order to 
investigate the sensitivity of our results to the choice of 
energy loss parameters in the model of [16], we consider 
three cases of parton energy loss. First we take energy 
loss to be constant in each parton scattering, = const., 
then we consider LPM energy-dependent energy loss [15], 
€ a n = a s \J E^Etpm and BH energy-dependent energy 
loss, e° = kE^. The scale separating the LPM region of 
validity from the BH region is given by the product of 
the average < p\ >, pr being the kick that parton gets 



from each collision and the mean free path, A Q . Taking 
< p T > to be lGeV 2 and A a to be 1 fm, one gets the 
scale to be about 5 GeV. Above this scale, one would ex- 
pect LPM effect to be dominant. However, increasing the 
value of A a or < p\ > shifts the transition place to higher 
values of px- Furthermore, one should keep in mind that 
LPM effect in QCD has been derived for static scatterers 
[15], which may not be suitable approximation in case on 
hot QGP. We find that the parton energy loss is mostly 
responsible for the observed suppression of hadron spec- 
tra in heavy ion collisions and that nuclear shadowing 
is a small effect (< 10% at most). Furthermore, the pr 
spectra are quite insensitive to the choice of scales [10]. 
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FIG. 1. Inclusive 7r° cross section in proton-proton and in 
Au-Au collisions at ^fs = 200 GeV. 

We show our results for the inclusive 7r° production in 
Au-Au collisions at ^/s = 200 GeV in Fig. 1 for different 
choice of parton energy loss parameter, e°. We chose 
the value for the constant energy loss, e° = 0.3GeV, 
Elpm — 0.35GeV, and the fraction of the energy lost 
per scattering, k — 6%, such that the suppression at 
Pt = 3 GeV is consistent with the observed suppression 
for pt < 3GeV [17]. Then we predict the shape of pt 
distribution for pt > 3 GeV for different choices of en- 
ergy loss. We note that pt distribution of neutral pions 
is much steeper in case of BH energy-dependent parton 
energy loss. 

In Fig. 2 we show the ratio of 7r° production in Au-Au 
collisions to the one in p-p collisions for different choice 
of parton energy loss parameter, e^. We find that for 
constant energy loss and for LPM energy-dependent en- 
ergy loss, the ratio increases with pt, while for the BH 
case, the ratio slightly decreases with pt- We show that 
the ratio is very sensitive to the fraction of energy loss, 
for example, for k = 10%, it decreases from about 50% 
to 35% for pt increasing from 3 GeV to 12GeV, while for 
k = 3% it decreases from 90% to 96%. When n = 6%, 
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the suppression increases from 76% at pt — 3 GeV to 
84% at p T = 8GeV, in agreement with recent PHENIX 
data [18]. 
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FIG. 2. Ratio of inclusive n cross sections in Au-Au colli- 
sion and in proton-proton collision at i/s = 200 GeV. 



production cross section reliably. Here we present results 
for prompt photon production in proton-proton and in 
Au-Au collisions at RHIC energy of 200GeV using the 
NLO pQCD code of Aurenche et al. with nuclear shad- 
owing and medium induced parton energy loss effects in- 
cluded. We use the same initial parton distributions, 
choice of scales, nuclear shadowing function and choice 
of parton energy loss, e° = 0.06i?° as for pions. However, 
we note that in prompt photon production there are two 
types of subprocesses that contribute, direct processes 
and bremsstrahlung processes, the later one being con- 
voluted with the fragmentation function [9] . Fragmenta- 
tion functions, D c /^(z,Q *), that describe the transition 
of the partons into the final-state 7, without medium 
effects, are extracted from e + e~ data [20]. In Fig. 3 
we show our results for prompt photon distribution in 
proton-proton and in Au-Au collisions at y/s = 200GeV. 
We note that suppression of prompt photons produced 
in heavy- ion collisions is much less that for 7r°. This is 
due to the fact that only bremsstrahlung processes are 
affected by the parton energy loss, which contribute 24% 
to the cross section at p T — 3 GeV and 6% at pt = 12 
GeV. For the same reason, prompt photons are not very 
sensitive to a different choice of parton energy loss. 



Clearly measurements of inclusive pion production at 
RHIC energies for pt > 3 GeV have provided valuable in- 
formation about the medium induced parton energy loss 
since nuclear shadowing effects are very small (few %) 
and the observed suppression of 7r° in Au-Au collisions 
can be described with an energy-dependent parton en- 
ergy loss. However, theoretical expectations that a LPM 
type of energy loss should be present is clearly not seen 
for pt < 5 — 6 GeV while LPM type energy loss seem 
to work for higher pr- One should keep in mind that 
all theoretical calculations of LPM energy loss are for 
very energetic partons and at high pr and applicability 
to RHIC at low to intermediate pr's is not appropriate. 
BH type of energy loss seem to reproduce the data quite 
well up to pt ~ 8 — 10 GeV but will fail beyond that. 
The most interesting aspect of this is that there seems 
to be different energy loss mechanisms at work at differ- 
ent pr's. Understanding ir° production at RHIC is also 
important because large-p-r pions form a significant 
background for the prompt photons. 

Studying photon production at the RHIC is also of spe- 
cial interest, as it has been suggested as an elegant sig- 
nal for detecting the formation of a quark-gluon plasma 
(QGP) in heavy-ion collisions [19]. Photons can be pro- 
duced at different stages of the heavy-ion collision: pho- 
tons produced at the early stages of the collision through 
hard parton-parton scatterings or they can be emitted 
from a thermalized quark gluon plasma or hadronic gas. 
Furthermore, prompt photons are an important back- 
ground to thermal photons in the low to intermediate pt 
region and are dominant in the high p t region. Therefore, 
it is extremely important to be able to calculate their 
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FIG. 3. Prompt photon production in proton-proton and 
in Au-Au collisions at y/s = 200 GeV. 



In Fig. 4 we show the ratio of prompt photon produc- 
tion in Au-Au and in proton-proton collision compared 
to the same ratio in case of 7r° production. We take e° to 
be energy-dependent, e° = 0.06£"°. We find that the ra- 
tio for prompt photon production is not very sensitive to 
the choice of energy loss, and for all choices considered, 
constant parton energy loss and energy-dependent par- 
ton energy loss, the suppression decreases with increasing 
Pt- We find that while the ir° suppression increases from 
76% to 88% for p T increasing from 3 GeV to 12 GeV, 
prompt photon suppression decreases from 40% to 15% 
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in the same range of px- 
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FIG. 4. Ratio of prompt photons produced in Au-Au col- 
lisions and in proton-proton collision at yfs = 200 GeV, com- 
pared with the same ratio for neutral pions. 

In Fig. 5, we show the ratio of prompt photons to 
7T° as a function of transverse momentum for ^fs = 200 
GeV, including nuclear shadowing and energy-dependent 
parton energy loss, = 0.06£^. We also show this ratio 
for proton-proton collisions at the same energy. In Au-Au 
collisions at RHIC, because of the large 7r° suppression 
relative to prompt photons at RHIC energies, this ratio 
increases with px approaching 1 at px ~ 10 GeV. 
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FIG. 5. Ratio of prompt photon and ty° cross sections at 
y/s = 200 GeV. 

In summary we have investigated inclusive pion and 
prompt photon production in proton-proton and in Au- 
Au collisions at RHIC. We have used the NLO pQCD 
results of Aurenche et al. for pion and prompt photon 
production in proton-proton collisions and included nu- 



clear shadowing and medium induced parton energy loss 
by modifying the final state pion and photon fragmenta- 
tion functions. 

We find nuclear shadowing effects to be rather small 
while medium induced parton energy loss results in large 
suppression of ir° production in Au-Au collisions at ^fs — 
200 GeV. We show that the suppression is sensitive to 
the type of energy loss assumed. In the case of constant 
or LPM energy loss, the suppression decreases with in- 
creasing px- On the other hand, when partons travers- 
ing the medium lose constant fraction of their energy, 
the suppression becomes stronger at higher px, which 
seems to be in agreement with recent PHENIX data [1]. 
The fact that one can reproduce the data with our sim- 
ple and quite naive model of energy loss is remarkable. 
Indeed, there are now more sophisticated models which 
include effects such as intrinsic transverse momentum of 
incoming partons, px broadening due to multiple parton 
scatterings as well as effects due to finite energy of the 
partons. These models also lead to the suppression of 
pion production in Au-Au collisions, however, the sup- 
pression is either decreasing with px [21] or is almost px 
independent for 2GeV < px < 20GeV at RHIC energies 
[22]. The effects of dynamical expansion of the collision 
region have been recently studied and related via scaling 
law to an equivalent static scenario [23]. The advantage 
of our simplistic model is that it can reproduce the data 
with only a few model-dependent assumptions. 

Using our results for pion production, we predict the 
prompt photon production cross section in proton-proton 
and Au-Au collisions at RHIC, where we use an energy- 
dependent parton energy loss, = 0.06i?^ and EKS 
nuclear shadowing function. We also present the ratio of 
prompt photon and 7r° production at RHIC, of relevance 
to separating different sources of photon production. At 
RHIC energy we find that at px — 3GeV direct processes 
contribute about 75%, while bremsstrahlung processes 
are remaining 25%. This changes significantly as one goes 
to higher energies, such as LHC, where bremsstrahlung 
contribution becomes dominant and photons become as 
suppressed as pions [10]. 
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